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Resumen 


Los datos obtenidos mediante satélites acerca de la precipitación pueden 
ser útiles en regiones donde ésta es muy heterogénea, como, por 
ejemplo, en los Andes, pues por lo general son escasamente 
monitoreadas. El objetivo de este estudio es precisar las características 
de las principales variables hidrológicas, y comprender el potencial de 
datos de precipitación estimados a partir del satélite en modelación 
hidrológica. Este artículo evalúa datos de precipitación obtenidos a partir 
de pluviómetros y satélites (productos TMPA V7, TMPA RT del satélite 
TRMM), como entrada al modelo hidrológico agregado GR2M, para 
simulación de caudales mensuales, en el periodo 2011-2015, en la 
cuenca del río Ilave, en el Altiplano peruano. Los resultados mostraron 
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un claro déficit hídrico de caudales debido a la disminución de las 
precipitaciones en el periodo de avenidas (“ 50%); además de un 
aumento de la evapotranspiración en el periodo de estiaje (“ 24%). Los 
resultados del análisis de datos de precipitación también indicaron que 
los TMPA V7 presentan mayor similitud que TMPA RT, con respecto a la 
precipitación observada en el periodo de avenidas. Del mismo modo, los 
resultados de modelación hidrológica con GR2M indicaron un mayor 
rendimiento a partir de la lluvia observada (pluviómetros) y TMPA V7, 
los cuales presentan una mayor similitud con caudales observados; a 
diferencia de los TMPA RT, es posible que los inconvenientes en el 
rendimiento se debieran a una inadecuada lluvia en el balance hídrico. 


Palabras clave: precipitación, modelación hidrológica, Andes, TRMM, 
satélite. 


Recibido: 23/05/2017 
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Introducción 


En las últimas décadas, numerosos datos globales y cuasiglobales de 
precipitación basada en satélites, observaciones terrestres o empleados 
de salidas de modelos de circulación general se han desarrollado en 
diferente resolución temporal (Huffman et a/., 2007). La precipitación es 
un componente importante del ciclo del agua y una entrada clave para 
las aplicaciones de modelado hidrológico. La cuantificación fiable de la 
distribución espacio-temporal de la lluvia es fundamental para el análisis 
de los fenómenos hidrológicos extremos, como sequías o inundaciones, 
en tiempo real o casi real. No obstante, las redes de medición 
(meteorológicas e hidrológicas) en regiones del Altiplano pueden ser 
escasas o incluso inexistentes. Además, la existencia de los Andes 
contribuye a la alta variabilidad espacial de las precipitaciones (Espinoza 
et al., 2009), manifestándose en que grandes cantidades de lluvia 
pueden precipitarse en pocos días, lo que incrementa el riesgo de 
inestabilidad del suelo, deslizamientos o inundaciones (Coscarelli 8 
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Caloiero, 2012; Zubieta € Saavedra, 2013; Zubieta, Saavedra, Silva, 8 
Giraldez, 2017a). Los productos de precipitación basados en satélites 
han facilitado cada vez más la estimación de precipitaciones en gran 
parte del mundo y se han convertido en una alternativa complementaria 
para aplicaciones hidrometeorológicas y estudios sobre el clima en 
regiones no monitoreadas. Entre ellas se encuentran los datos del 
Tropical Rainfall Measuring Mission (TRMM), Análisis Multi-Satélite de 
Precipitación (TMPA, por sus siglas en inglés) de Datos (Huffman, Adler, 
Bolvin, € Nelkin, 2010). TRMM fue una misión que comenzó en 1997, 
entre el National Space Development Agency del Japón (JAXA) y la 
National Aeronautics and Space Administration (NASA) de los Estados 
Unidos de América, la cual llegó a su fin el 8 de abril de 2015, como 
parte de su desactivación, al haber agotado sus reservas de energía 
(https://pmm.nasa.gov/trmm/mission-end). De forma reciente, la 
utilidad de estos datos en cuencas andino-amazónicas ha sido evaluada 
mediante la comparación de datos estimados y observados tanto de 
lluvia como caudal, dando a conocer que datos de precipitación a partir 
de satélites pueden ser considerados como entrada, para la modelación 
de caudales en cuencas (Zulkafli et a/., 2014; Zubieta, Getirana, 
Espinoza, €: Lavado, 2015; Milewski, Elikadri, £€ Durham, 2015; Zubieta, 
Getirana, Espinoza, Lavado-Casimiro, 8 Aragon, 2017b). Zulkafli et al. 
(2014) documentaron que el rendimiento obtenido de la modelación 
hidrológica tiende a ser insatisfactorio sobre regiones andino- 
amazónicas del norte del Perú y el sur del Ecuador. No obstante, se 
obtienen mejores resultados sobre regiones andino-amazónicas del sur 
del Perú, los cuales están asociados con estimaciones adecuadas de 
lluvia a partir de satélite (Zubieta et a/., 2015). 


Para evaluar la utilidad de los datos de precipitación basados en 
satélites es necesario investigar heterogeneidades temporales de la 
lluvia y sus impactos en el ciclo hidrológico. La variabilidad de los 
caudales en cuencas andinas es consistente con la variabilidad de las 
precipitaciones a escala estacional e interanual. Esto se presenta sobre 
todo en los Andes peruanos, que se caracteriza por diferentes regímenes 
de lluvias, debido a su extensión latitudinal y la presencia de la 
cordillera (Espinoza et a/., 2009). Por ejemplo, las precipitaciones en la 
cuenca del río Ramis en el Altiplano peruano presentan tendencias 
negativas de 0.70 mm/año, mientras las temperaturas máximas 
presentan tendencias positivas significativas promedio de 0.04 *C/año; 
asimismo, las temperaturas mínimas a nivel de cuenca muestran un 
incremento de 0.0004 “C/año (Belizario, 2015). Por otro lado, en los 
últimos años, El Altiplano ha sufrido graves eventos hidrológicos 
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extremos, como intensas sequías en 2008 y 2010 (ANA, 2010), e 
inundaciones en 2012 (RPP Noticias, 2012). 


El objetivo principal de este estudio es evaluar la capacidad de datos 
TRMM, productos TMPA V7 y TMPA RT (real time) de precipitaciones 
basadas en satélites, para representar campos de precipitación de 
utilidad, para la modelación hidrológica de la cuenca del río Ilave (CRI). 
La CRI se ubica en los Andes, del sur del Perú, en la región norte del 
Altiplano (69.4 oeste-70.38 oeste; 17.14 sur-15.8 sur), con área de 
drenaje de 7 790 km?*; con altitudes que van de 3 805 a 5 400 m de 
altitud, y una altitud media de 4 300 msnm (Figura la-b). La 
precipitación anual media a lo largo de la CRI es de entre 494 y 888 
mm/año, además presenta un caudal promedio de 32 m*/s. La amplitud 
del ciclo anual de precipitación es relativamente grande, con valores 
máximos entre diciembre y febrero, y mínimos entre junio y agosto. 
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Figura 1. a) Ubicación de estaciones climáticas; b) ubicación de la cuenca del río 
Ilave. 


También se consideraron los datos de precipitación observada, para 
efectos de comparación con las estimaciones satelitales. Cada conjunto 
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de datos de precipitación (observado y satelital) se utilizan como 
entrada para el modelo GR2M (Niel, Paturel, € Servat, 2003), el cual se 
basa en el modelo GR2M (Edijatno €: Michel, 1989; Kabouya, 1990), que 
fue revisado por Makhlouf y Michel (1994), para simulaciones de flujo 
fluvial y su comparación con caudales observados. GR2M se ha 
empleado para simular el comportamiento hidrológico de numerosas 
cuencas andino-amazónicas en el Perú (Lavado, Labat, Guyot, Ronchail, 
8 Ordoñez 2009; Lavado, Labat, Guyot, € Ardoin-Bardin, 2011; Vera € 
Felipe, 2011). Recientemente, Lujano, Felipe, Lujano y Quispe (2015) 
demostraron que los datos TMPA son factibles de emplear como entrada 
para el modelo GR2M en la cuenca del río Ramis, la cual se ubica al 
noroeste de la CRI. 


Datos 


Se recopilaron datos de precipitación diaria de 11 estaciones de lluvia 
del Servicio Nacional de Meteorología e Hidrología (SENAMHI). El 
periodo de datos seleccionado para el análisis hidrológico fue entre 2010 
y 2015. No obstante, para garantizar la máxima disponibilidad de datos 
pluviométricos e hidrométrica se seleccionó un periodo común (2011- 
2015) (Figura 1b, Tabla 1). Para evaluar su calidad, esta base de datos, 
compuesta por un total de 10 estaciones de lluvia sobre una base 
mensual, se sometió al método de vector regional (RVM, por sus siglas 
en inglés), que utiliza el concepto de precipitación media extendida al 
periodo de estudio (Hiez, 1977; Brunet-Moret, 1979; Vauchel, 2006). 


Tabla 1. Características generales de estaciones climáticas. 
Latitud, longitud, altitud, lluvia, temperatura máxima (T máx), temperatura mínima (T 
mín) y porcentaje de la disponibilidad de datos en el periodo 2011-2015. 


Altitud | Lluvia T T Rara 

Nombre Latitud Longitud | (msnm) (mm) máx mín datos (%) 
Capazo -17.19 -69.74 4 530 464 13.6 -6.3 99.8 
Mazo Cruz -16.74 -69.72 4 003 564 17.5 -5.0 99.9 
Laraqueri -16.15 -70.07 3 900 788 16.7 -1.6 99.8 
Los Uros -15.82 -69.94 3 808 810 15.8 3.9 99.9 
89 
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Rincón La Cruz -15,99 -69.81 3 935 872 15.3 2.7 99.9 
llave -16.07 -69.66 3871 647 15.3 1.5 99.9 
Callacame -16.57 -69.32 3951 413 14.0 1.0 99.8 
Pizacoma -16.91 -69,37 3 930 583 17.0 -0.9 99.9 
Ichuña -16.13 -70.55 3800 595 20.1 2.5 99.9 
Calacoa -16.74 -70.68 3 260 386 25.0 4 99.8 


El método de mínimos cuadrados se aplicó para encontrar un índice 
pluviométrico anual regional Zí y una precipitación media extendida. 
Este cálculo puede realizarse minimizando la suma de la ecuación (1), 
en donde P,; representa la lluvia anual para la estación j, P, es la 
precipitación media extendida; ¡, el índice años; j, el índice de la 
estación; N, el número de años, y M es el número de estaciones. Por 
último, la serie de datos de Zí se llama vector regional anual de los 
índices pluviométricos: 


Y) (2-2) (1) 


De este modo, se consideró una misma zona climática de la CRI, 
experimentando un mismo régimen de lluvias; se asumió que la 
precipitación anual en las estaciones de la cuenca mostró la 
proporcionalidad entre las estaciones; TRMM provee datos derivados 
entre sus productos, los cuales son evaluados en este estudio tanto: 


1. TMPA 3B42 versión 7, que se obtiene a partir de un conjunto de datos 
de muchos sensores a bordo de satélites de precipitación tanto en 
tiempo real y casi-real (datos TMPA 3B42 están disponibles en 
https://pmm.nasa.gov/data-access/downloads/trmm). 


TMPA muestra la experiencia más exitosa, ya que su procesamiento ha 
sido perfeccionado, al emplear datos de mejor escala en relación con 
otras fuentes de datos, siendo primordial crear promedios a nivel 
espacio-temporal que sean apropiados para la aplicación del usuario 
(Huffman et a/., 2010). 


2. Otro conjunto de datos son las estimaciones de precipitación TMPA- 
RT. Estos datos en relación con TMPA V7 no incorporan mediciones de 
calibración de estaciones de lluvia, las cuales son incorporadas más de 
un mes después que de los datos de satélite. 
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Tanto TMPA V7 y TMPA RT se caracterizan por proporcionar datos de 
precipitación con resolución espacial de 0.25% * 0,25% y resolución 
temporal de tres horas. 


Se recopilaron datos de precipitación diaria de 11 estaciones de lluvia 
del Servicio Nacional de Meteorología e Hidrología (SENAMHI). El 
periodo de datos seleccionado para el análisis hidrológico fue entre los 
años 2010 y 2015. No obstante, para garantizar la máxima 
disponibilidad de datos pluviométricos e hidrométrica, se seleccionó un 
periodo común (2011-2015) (Figura 1b, Tabla 1). 


Para evaluar su calidad, esta base de datos compuesta por un total de 
10 estaciones de lluvia sobre una base mensual, la cual se sometió al 
método de vector regional (RVM), que utiliza el concepto de 
precipitación media extendida al periodo de estudio (Hiez, 1977; Brunet- 
Moret, 1979; Vauchel, 2006). El método de mínimos cuadrados se aplicó 
para encontrar un índice pluviométrico anual regional Zi, y una 
precipitación media extendida. Para este fin, el cálculo puede realizarse 
minimizando la suma de la ecuación (1), en donde P,; es la precipitación 
anual en la estación j y P, es la precipitación media extendida; ¡, el 
índice años; j, el índice de la estación; N, el número de años, y M es el 
número de estaciones. Por último, la serie de datos de Zi se llama 
vector regional anual de los índices pluviométricos. 


Metodología 


El GR2M es un modelo agregado que funciona a paso de tiempo mensual 
y requiere información promedio de variables climáticas para toda el 
área de drenaje hasta la estación Puente Carretera llave (Figura la). 
Ante el número reducido de estaciones (10) e incertidumbre en algunas 
regiones no monitoreadas de la cuenca, se procedió a interpolar datos 
de precipitación a la misma resolución de TMPA (0.25% * 0,259) 
mediante el método Kriging, aprovechando la interrelación que hay 
entre las estaciones ubicadas dentro y fuera de la cuenca. El método 
seleccionado en este estudio es el Kriging ordinario, que consiste en 
cuantificar el supuesto de que las mediciones cercanas tienden a ser 
más similares que otras más alejadas. La información detallada acerca 


91 
Tecnología y ciencias del agua, 9(5), 85-105, DO!:10.24850/j-tyca-2018-05-04 


Tecnología y e 


CienciaszAgua 


del método de interpolación Kriging puede encontrarse en Lichtenstern 
(2013). Finalmente, se determinaron los valores promedio de 
precipitación para cada mes y para toda la cuenca durante el periodo 
2010-2015. 


Con base en la disponibilidad de datos de temperatura máxima y 
mínima, la variable evapotranspiración (ETP) se estimó con el modelo de 
Hargreaves y Samani (1985), cuyo cálculo se muestra en la ecuación 


(2): 
ETP = 0.0023 (Tin + 17.8 Tmax — Tmin) > * Ra (2) 


donde ETP es evapotranspiración (mm/día); T,,, temperatura media 
(*C); Tmaxr, temperatura máxima (*C); Tin, temperatura mínima (*C), y 
Ra es radiación extraterrestre (mm/día). Este método proporciona una 
apropiada aproximación con respecto a datos obtenidos por lisimetros 
en regiones del norte del Altiplano (mayor humedad), siendo éstos 
similares a aquellos obtenidos por la ecuación de FAO Penman-Monteith 
(García, Raes, Allen, 8. Herbas, 2004). De hecho, estudios realizados en 
el Altiplano boliviano empleando variedades de pastizal (por encima de 3 
000 msnm) como cultivo de referencia sugieren que los métodos de 
Penman-Monteith y Hargreaves-Samani concuerdan de modo apropiado 
con datos a partir de lisímetros (García, Raes, € Jacobsen, 2003; García 
et al., 2004; Vacher, Imana, € Canqui, 1994). No obstante, Hargreaves 
y Samani (1985) tienden a ligeramente subestimar datos observados en 
la región sur del Altiplano, donde los factores aerodinámicos que afectan 
la evapotranspiración son más importantes que en el norte del Altiplano 
debido a la llanura y a un aumento de la aridez (García et al., 2004). 


GR2M depende de dos parámetros: X1 (máxima capacidad de 
almacenamiento del reservorio en mm) y X2 (coeficiente de intercambio 
de agua subterránea, adimensional). Conceptualmente, el modelo se 
basa en la transformación de la precipitación-escorrentía mediante dos 
funciones: de producción y de transferencia. El modelo GR2M es un 
modelo a dos reservorios, donde la función de producción se organiza 
alrededor de un reservorio llamado reservorio-suelo; la función de 
transferencia es regida por el segundo reservorio, llamado reservorio de 
agua gravitacional, donde el aporte es instantáneo al inicio del paso de 
tiempo, luego el reservorio se vacía de manera gradual. El nivel de este 
reservorio determina el caudal que puede liberar (Figura 2). 
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Si S+ XP with p= cnt] 3 
ETP P beats Xx 3) 
Vo (3) a (4) 
(4) P=P+8S-S, 
Evapotranspiración 
Ñ S,= S¡(1-y) with ve (5) 
Ñ S X; 
Reservorio de 1+ ob = z 
Producción x Ss X; 
S, P, =S,-S (6) 


A 3113 
(e) P, P; 12) 
—— X 
P: 2 (7) 


(7) 8 P;, = P, + P, 
| de (8) 
Intercambio Xx PP R¡ =R+H 
Subterráneo R 
(9) Reservorio R>X.R; (9) 
de agua R? 
(10) 10) gravitacional. O= Fr R=R,-Q (10) 


Figura 2. Modelo GR2M de Mouelhi (2003). 
Fuente: Cemagref (www.cemagref.fr). 


Las variables de entrada (P) y salida (ETP) del reservorio-suelo 
(capacidad máxima X1) se calculan con base en el almacenamiento del 
reservorio de producción (S). Una parte de la precipitación P1 hace 
referencia a la diferencia entre P y la precipitación que ingresa en el 
suelo, y otra parte P2 hace referencia a la cantidad de agua que 
contribuye al intercambio subterráneo. La suma tanto de Pi y P2 
corresponde a la precipitación efectiva P3, que ingresa en el reservorio 
de agua gravitacional (R) (capacidad máxima fijada en 60 mm). El 
intercambio subterráneo (R2) se incluye en el modelo X2. Por último, 
como resultado para el caudal de salida, el reservorio de agua 
gravitacional se libera de acuerdo con una función cuadrática. La 
optimización del modelo responde a dos parámetros: X1, capacidad del 
reservorio-suelo en milímetros; y X2, coeficiente de intercambios 
subterráneos (adimensional). Para ello, se realizó una calibración 
manual previa, para emplear la herramienta SOLVER de Excel. 


Para evaluar el rendimiento del modelo se utilizaron las siguientes 
funciones objetivo: coeficiente de Nash-Sutcliffe (NS); diferencia entre 
volúmenes calculados y observados mediante el error en el volumen 
(AV), y coeficiente de determinación mostrados (po,,..o..,) en las 


ecuaciones (11), (12) y (13): 
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21 Qops(t) => Oops)? 
A= »(Qops(t)) E »(Qops(t)) (12) 
N(Qops(t)) 
Pdobs ¡Qcal = (13) 


Donde Qo»s es el caudal observado; Q«a1, el caudal estimado; 0... 001, la 
covarianza entre Qo»s Y Qca; por Último, 02,,. Y o,. representan la 
desviación estándar de las variables Qo»s Y Qca, respectivamente. 


Resultados y discusión 


En términos promedio, el análisis visual del régimen de 
evapotranspiración estimada para el periodo 2011-2015 no muestra 
variaciones significativas durante el periodo de avenidas (“ 166 mm, 
diciembre-febrero); no obstante, hay un incremento aproximado de 64 a 
80 mm durante el periodo de estiaje. Basado en el balance hídrico, esto 
se asocia con la reducción de la precipitaciones identificada a lo largo del 
periodo de estudio (“ 50%) (Figura 3a). La precipitación media anual en 
la CRI indica que TMPA V7 (836.4 mm) y TMPA RT (1 022.5 mm), 
sobreestiman datos observados (688.0 mm) en 21.6 y 48.6%, 
respectivamente. 
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AAN 


mu Precipitación (mm) 
—— ETP (mm) 


164.2 7.9 75.1 
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Figura 3. a) Comparación entre evapotranspiración y lluvia observada, diagramas de 
dispersión de los productos de lluvia; b) TMPA V7, y c) TMPA RT en relación con la 
lluvia observada; d) lluvia promedio mensual para cada producto de lluvia. 


El análisis de la precipitación mensual de la CRI muestra que las 
estimaciones de datos TMPA V7 son más realistas en relación con los 
datos observados al ser identificado un alto coeficiente de determinación 
entre ellos (1? = 0.86) y menor error estándar (EME = 25 mm) (Figura 
3b). Sin embargo, TMPA V7 muestra diferencias con respecto a datos 
observados en las estimaciones de lluvia por encima de 50 mm/mes. 
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Asimismo, datos TMPA RT muestra mayores diferencias que TMPA V7 
para estimaciones por encima de 50 mm/mes, a pesar de presentar un 
coeficiente de determinación relativamente aceptable (1? = 0.76) y un 
mayor error (EME = 33.1 mm) (Figura 3c). 


Para un análisis temporal de la lluvia se calculó la lluvia promedio 
mensual para cada conjunto de datos (pluviómetros, TMPA V7, TMPA 
RT). El análisis corroboró que las estimaciones TMPA V7 son las más 
aproximadas en relación con los datos observados durante el periodo de 
avenidas (diciembre-marzo) (Figura 3d). Pero TMPA V?7 muestra 
tendencia a sobreestimar la lluvia observada durante el periodo de 
estiaje en “ 22 mm (mayo-septiembre). Por otro lado, datos TMPA RT 
muestran tendencia a sobreestimar datos observados 
predominantemente a lo largo de todo el año. 


Las diferencias identificadas en las estimaciones satelitales con respecto 
a datos de lluvia observada pueden estar asociadas con variaciones 
observacionales y/o provenientes de datos de sensores en la entrada de 
datos al algoritmo utilizado para la estimación de precipitación satelital 
(Huffman et al., 2007). De hecho, se encontraron resultados similares 
en la cuenca del río Ramis, en el Altiplano peruano (Lujano et a/., 2015), 
u otras regiones de los Andes del Perú y Ecuador (Zulkafli et a/., 2014; 
Zubieta et al., 2015). 


Para el proceso de calibración del modelo se realizó una evaluación 
previa en la selección del número de años de los periodos de calibración 
y validación. En condiciones promedio, para un periodo de calibración de 
un año (periodo de validación cuatro años), el rendimiento del modelo 
tiende a ser mayor que al seleccionar 2 o 3 años (periodo de validación 
3 y 2 años), los cuales tienen reducciones de 30 y 45% en el 
rendimiento, respectivamente. Por tal motivo, para la calibración del 
modelo, tanto de lluvia observada y de satélite, se consideró sólo el 
periodo 2011, mientras que para la validación se consideró el periodo de 
2012-2015. 


Las variables de entrada al modelo fueron los valores promedio 
mensuales de precipitación, evapotranspiración potencial y caudales. El 
proceso de optimización se basó en la obtención de parámetros idóneos 
(X1 y X2), del modelo de acuerdo con la ecuación (10), empleando 
como variable de entrada al modelo, tanto datos de lluvia observada 
como datos satelitales (TMPA V7, TMPA RT). Los valores X1 y X2 
optimizados para cada conjunto de datos de precipitación se muestran 
en la Tabla 2, donde también se observan los valores encontrados para 
el coeficiente de Nash, desbalance hídrico y coeficientes de correlación, 
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al realizar la evaluación de datos de caudales modelados con respecto a 
caudales observados. 


Tabla 2. Parámetros de calibración, coeficientes de rendimiento Nash-Sutcliffe, 
diferencias entre volúmenes, y coeficiente de determinación entre caudales calculados 
y observados. 


Lluvia 
Parámetro observada |TMPA V7 |  TMPA RT 

X1: reservorio de producción (mm) 5.76 5.57 6.35 
X2: parámetro de intercambio (mm) 1.00 0.96 0.74 
Coeficiente de Nash-Sutcliffe (NS) 0.95 0.74 0.52 
% de diferencia entre volúmenes 

calculados y observados (AV) 5 23 41 
Coeficiente de determinación (R?) 0.85 0.64 0.45 


Similar al régimen de precipitaciones, el régimen de caudales 
observados muestra también un claro retroceso de 2011 a 2015 tanto 
en el periodo de avenidas o estiaje, y del patrón similar para cada serie 
de caudales estimadas a partir de datos de satélite (TMPA V7 y TMPA 
RT). De acuerdo con el balance hídrico y el alto rendimiento del modelo 
GR2M en la estimación de caudales empleando datos observados, este 
déficit se debe a la disminución predominante de las precipitaciones, el 
cual también es identificado en la serie de caudales obtenidos a partir de 
satélite. Del mismo modo, está influenciado por el incremento de la 
evapotranspiración estimada en estiaje durante el periodo 2011-2015 
(Figuras 4a y 3a). 
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Q(m?/s) P(mm) 


Ene 11 Ene 12 Ene 13 Ene 14 Ene 15 


150 | === Modelado 
250 | ——— Observado 


| Precipitación (mm) 


Ene 11 j Ene 12 Ene 13 j Ene 14 Ene 15 


Ene 11 Ene 12 Ene 13 Ene 14 Ene 15 


Figura 4. Caudales observados y modelados en la estación Puente Carretera llave 
desde enero de 2011 hasta junio de 2015. Los caudales modelados son calculados 
empleando datos de precipitación obtenidos a partir de a) pluviómetros; b) TMPA V7, y 
Cc) TMPA RT en la cuenca del río Ilave. 


La disminución de las precipitaciones puede estar asociada con el 
incremento sustancial en la ocurrencia de frentes cálidos en el Altiplano, 
los cuales pueden aumentar la evaporación atmosférica, además de 
disminuir la vegetación (Thibeault, Seth, 8 Garcia, 2010). 


El análisis visual de los hidrogramas revela que los caudales obtenidos a 
partir de lluvia observada (pluviómetros) muestran mayor similitud con 
caudales observados registrados en la estación de aforo Puente 
Carretera llave (NS = 0.95 y AV = -11.5 %), que se muestran en la 
Figura 4a. El régimen de precipitaciones es consistente en relación con 
los caudales obtenidos en el balance hídrico (ver Figura 4a). Las 
fluctuaciones estacionales en el régimen de caudales están muy bien 
representadas por el modelo GR2M, al emplear datos observados, 
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principalmente durante el periodo de validación, tal es el caso de los 
caudales registrados en el periodo de avenidas 2012, 2014 y 2015. 


No obstante, el rendimiento del modelo GR2M al utilizar datos TMPA V7 
en la CRI es mínimamente aceptable, NS = 0.74, AV = -7.3% a pesar 
de presentar un NS relativamente alto (Figura 4b), pues no tiene la 
capacidad de representar de manera adecuada todos los caudales picos, 
lo que podría deberse a la sobreestimación de la lluvia durante la 
temporada de inicio de máximas avenidas (diciembre-enero). 


Sin embargo, puede reproducir con mejor aproximación caudales 
estacionales, sobre todo durante los periodos de estiaje. Se puede 
observar menor rendimiento al obtenido con TMPA V7 cuando se utiliza 
TMPA RT (NS = 0.52, AV = -28%) (Figura 4c). Basado en los resultados 
adecuados en la modelación, a partir de datos de lluvia observada, el 
muy bajo rendimiento al emplear datos TMPA RT se debe a una lluvia 
estimada por satélite inadecuada, caracterizada en valores promedio por 
una marcada sobrestimación con respecto a datos observados 
(pluviómetro) (31%) y datos TMPA V7 (35%) durante el periodo de 
avenidas (Figura 3d). Asimismo, esta sobrestimación también se aprecia 
durante el estiaje, la cual asciende hasta en 81% con respecto a la lluvia 
observada (pluviómetros) y 16% con respecto a la lluvia TMPA V7. 


Los caudales mensuales 2011-2015 obtenidos a partir del satélite se 
compararon con datos observados provenientes de la estación Puente 
Carretera llave mediante diagramas de dispersión (Figura 5a-c), y 
muestran que las series generadas por el modelo que usa datos 
observados (pluviómetros) y TMPA V7 representan mejor los caudales 
observados, al presentar coeficientes de determinación aceptables (0.85 
y 0.64, respectivamente). Sin embargo, los caudales (Q > 75 m?/s) 
presentes de manera predominante durante el periodo de máximas 
avenidas (enero-marzo) presentan una mayor dispersión de los datos, a 
diferencia de otros periodos (estiaje durante Junio-Julio o inicio de 
precipitación entre Setiembre-diciembre). En general, esto sugiere una 
mejor aproximación de datos TMPA V?7, siendo éstos los de mayor 
similitud a caudales observados. No obstante, los caudales simulados 
usando TMPA RT muestran un bajo coeficiente de determinación (1? = 
0.45). 


99 
Tecnología y ciencias del agua, 9(5), 85-105, DO!:10.24850/j-tyca-2018-05-04 


Tecnología y 


CienciaszAgua 
PLU = 0.89'0o - 1.21 250 RT =0.52*00 + 7.69 
a o 
E 200 
a 
S 150 
3 100 
> | 
S 
EME=18.8 di EME=36.5 
r? = 0.85 0 ld r? = 0.45 
O 50 100 150 200 250 O 50 100 150 200 250 O 50 100 150 200 250 
Caudal observado (Qo) Caudal observado (Qo) Caudal observado (Qo0) 


Figura 5. Caudales observados y modelados empleando a) Lluvia observada y lluvia 
(PLU) obtenida por satélites b) TMPA V7, c) TMPA RT. 


La mayor similitud entre datos TMPA V7 y datos observados de lluvia, 
además de la sobrestimación de TMPA RT con respecto a datos 
observados (Figura 4b-c), son congruentes con el impacto reflejado en 
el parámetro de máxima capacidad de almacenamiento del reservorio 
X1, pues estos parámetros obtenidos en el proceso de optimización 
muestran ser más similares cuando son provistos por TMPA V7 y datos 
observados (Tabla 2). 


Por otro lado, los caudales modelados con GR2M empleando lluvia 
satelital (TMPA V7, TMPA RT) tienden a subestimar caudales 
observados. Esto se asocia principalmente con el proceso de 
optimización del modelo, pues un menor intercambio subterráneo (X2) 
es efectivo cuando se modela con lluvia satelital provista por TMPA V7 
OQ = 0.96) y TMPA RT (X2 = 0.74) (ver Tabla 2). Esto sugiere un 
impacto negativo en la producción final de caudales a medida que X2 
sea menor, lo cual implica que aunque la precipitación satelital tienda a 
sobrestimar lluvia observada en el periodo de avenidas, los resultados 
del modelo empleando esta lluvia tienden a subestimar los caudales 
observados. 


Conclusiones 
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Se ha identificado un claro déficit en el régimen hídrico de caudales del 
río Ilave en el Altiplano peruano para el periodo 2011-2015. De acuerdo 
con el balance hidrológico realizado, el déficit hídrico se debe a la 
disminución de las precipitaciones en periodo de avenidas (“ 50%) 
además de un aumento de la evapotranspiración en el periodo de estiaje 
( 24%). 


Dos conjuntos de datos de lluvia mensual basados en el satélite TRMM 
(TMPA V7, TMPA RT) se compararon con datos de precipitación 
observada, basada en pluviómetros sobre la cuenca del río Ilave. En 
condiciones promedio, y a partir de análisis estadístico, TMPAV?7 
presenta una mayor asociación con respecto a la precipitación 
observada, principalmente durante el periodo de avenidas (enero- 
marzo); no obstante, TMPA V7 muestra tendencia a sobreestimar lluvia 
observada durante el periodo de estiaje (mayo-setiembre). Por otro 
lado, datos basados en TMPA RT muestran de manera predominante la 
sobreestimación a lo largo del año. 


Para investigar la ventaja de datos satelitales en hidrología, los 
conjuntos de datos de precipitación observada basada en pluviómetros, 
TMPA V7 y TMPA RT se emplearon como variable de entrada al modelo 
hidrológico GR2M. El análisis de los hidrogramas revela que los caudales 
obtenidos a partir de lluvia observada (pluviómetros) muestra mayor 
similitud con caudales observados en la estación de aforo Puente 
Carretera llave (NS = 0.95). Cabe resaltar que el rendimiento del 
modelo GR2M al emplear datos TMPA V7 puede representar de forma 
adecuada el ciclo estacional de caudales (NS = 0.74), pero no presenta 
habilidad para representar de modo correcto caudales pico. Esto sugiere 
una inadecuada estimación de lluvia durante el mes de febrero, 
caracterizada por una subestimación de lluvia. 


Se encontró un menor rendimiento al obtenido con TMPA V7 al usar 
TMPA RT (NS = 0.54). Con base en resultados adecuados en la 
modelación, a partir de datos de lluvia observada, el muy bajo 
rendimiento al emplear datos TMPA RT se debe a una inadecuada lluvia 
estimada por satélite, caracterizada por una marcada sobrestimación 
con respecto a datos observados. De hecho, los caudales muestran que 
las series generadas por el modelo a partir de datos observados 
(pluviómetros) y TMPA V7 representan mejor los caudales observados 
durante el periodo de avenidas. 


Errores en las simulaciones de caudal de GR2M pueden estar asociados 
principalmente con la entrada de datos de lluvia o evapotranspiración, 
además de la representación limitada de procesos físicos del modelo. Sin 
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embargo, los resultados muestran que es posible utilizar datos 
obtenidos a partir de satélites y usarlos en modelación hidrológica de 
pequeñas cuencas en los Andes altiplánicos del Perú, para simular de 
modo apropiado los caudales. 


Evaluar los datos de la nueva generación de estimaciones de lluvia 
proporcionada por el satélite GPM (Global Precipitation Measurement) 
(Schwaller 8 Morris, 2011) es un área de estudio interesante para la 
hidrología. Este tipo de datos provee la posibilidad de tomar ventaja de 
datos subdiarios o subhorarios de precipitación para estimación de 
caudales de la misma resolución temporal usando modelos hidrológicos. 
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